Testing tripartite Mermin inequalities by spectral joint-measurements of qubits 
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It is well known that Bell inequality supporting the local realism can be violated in quantum mechanics. 
Numerous tests of such a violation have been demonstrated with bipartite entanglements. Using spectral joint- 
measurements of the qubits, we here propose a scheme to test the tripartite Mermin inequality (a three-qubit 
Bell-type inequality) with three qubits dispersively-coupled to a driven cavity. First, we show how to generate 
a three-qubit Greenberger-Horne-Zeilinger (GHZ) state by only one-step quantum operation. Then, spectral 
joint-measurements are introduced to directly confirm such a tripartite entanglement. Assisted by a series of 
single-qubit operations, these measurements are further utilized to test the Mermin inequality. The feasibility of 
the proposal is robustly demonstrated by the present numerical experiments. 

PACS number(s): 03.65.Ud, 42.50.Dv, 42.50.Pq 
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I. INTRODUCTION 

Entanglement 01 is at the heart of the quantum theoryand 
the crucial resource of quantum information processing [0 . 
It is also one of the most important ingredients of various in- 
triguing phenomena, e.g., quantum teleportation HID], secret 
sharing [:6:], and remote state preparation [7], etc. Therefore, 
generating and verifying the existence of entanglements are of 
great importance. 

Since Bell inequality H] and its CHSH version J^] was for- 
mulated to test the correlations between two particles, nu- 
merous experiments with bipartite entanglement, e.g., pho- 
tons ifioll . tra ppe d ions ITTl H2I1 . neutrons lll3ll and Josep- 
son junctions lll4l[l5ll . etc., have been demonstrated to probe 
the nonlocal nature of quantum mechanics. As all these ex- 
periments support quantum mechanics and rule out the local 
hidden-variable theories, Bell inequality can be served as an 
important witness of quantum entanglement. 

With the developments of quantum technology, entangle- 
ment shared by multiple particles play more and more impor- 
tant roles for large-scale quantum information processing and 
many-body quantum mechanics. Experimentally, multipartite 
entangled states have been demonstrated with photons lTl6l — 



I19H . trapped ions 1 2C -12211. Rydberg atoms 023H . and also 
Josephson circuits 11241 12511 . etc. Basically, multipartite en- 
tanglement can be robustly verified by the standard quantum- 
state tomographic technique, i.e., reconstructing their den- 
sity matrixes by a series of quantum measurements. Instead, 
one can also verify entanglement by testing the violation of 
the multipartite Bell-type inequality, such as Mermin inequal- 
ity JM]: 
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with three-qubit systems. Indeed, the violation of this inequal- 
ity has been experimentally demonstrated with three-photon 
entanglement OEl. Above, {9 U 9 2 , 9 3 , 0'i, 9 2 , 9' 3 } are 
the set of controllable local-variables of the three independent 
particles, and the correlation function E(9\, 9 2 , 9 3 ) is the en- 
semble average over the measurement outcomes for the local 
settings: 9\,9 2 ,9 3 . 

As a possible experimental demonstration, in this paper we 
discuss how to perform the test of a tripartite Mermin in- 
equality with three qubits coupled to a driven cavity. Two 
main contributions in the present proposal are: (i) an one-step 
approach is proposed to generate the desired Greenberger- 
Horne-Zeilinger (GHZ) state rf27ll . and (ii) a spectral measure- 
ment method is introduced to implement the joint measure- 
ments of these three qubits. In principle, our proposal could be 
further generalized to the cases with more than three particles. 
The paper is organized as: In Sec. II, we briefly describe how 
to generate the desired tripartite GHZ entangled state of three 
qubits coupled dispersively to a driven cavity. Then, by intro- 
ducing a spectral joint-measurement method via detecting the 
photon transmission through the driven cavity, we propose a 
simple two-step method to confirm such a tripartite GHZ en- 
tanglement. In Sec. Ill, we propose how to encode various 
local variables into the prepared GHZ entanglement via per- 
forming suitable single-qubit operations, and implement the 
test of the desired Mermin inequality by the introduced joint- 
measurements. Discussions on the feasibility of our proposal 
are given in Sec. IV. 



II. GENERATION AND CONFIRMATION OF THE GHZ 
STATE OF QUBITS COUPLED TO A DRIVEN CAVITY-QED 
SYSTEM 

A. Preparation of tripartite GHZ state by only one-step 
quantum operation 
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We consider a driven cavity-qubit system, wherein three 
qubits without interbit interaction are respectively coupled to 
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a common driven cavity. In principle, such a cavity-qubit sys- 
tem can be described by the Tavis-Cummings Hamiltonian 
lElll (h = 1 throughout the paper) 

H T c = w r a t a + ^ [-^-a Zj + g j (a^a- j + aa +j )], (2) 
j=i,2,3 

where aft' and o± } are the ladder operators for the photon 
field and the jth qubit, respectively; oj r is the cavity frequency, 
cjj the jth qubit transition frequency, and gj the coupling 
strength between the jth qubit and the cavity. The driving 
of the cavity can be modeled by 

H d = e(t)(a*e- iu * t + ae*"-*), (3) 

where e(t) is the amplitude and ujd the frequency of the exter- 
nal drive. 

Following Ref. ll29tl . after a displacement transformation 
D(a) = cxp(aa^ — ct*a), the displaced Hamiltonian of the 
composite system reads 

H T = D\a)(H TG + H d )D(a)-iD\a)D(a), 
= cj r a^a+ ^ l~^ az J +5j( at < T - J + arT +j) 

.7=1,2,3 

+g j (a*a- j + aa +j )]. (4) 

Now, we let LUd = ujj, a = —iuj r a — is exp(-iujdt), and 
work in a rotating frame defined by U\ = exp[— it{ui r a^a + 
Ud<J z /2)], the effective Hamiltonian of the cavity- 

qubit system takes the form 

St = ^ [ttjVxj +9j(a^-i exp(-iSt) 

3 = 1.2,3 

+acr +J exp (iSt))], (5) 

with the qubit-drive detuning 5 = Ud — u r and the Rabi 
frequency IX, = egj/S. Changing to the orthogonal bases 
|±j) = (|lj) ± \0j))/^/2, and in the interaction picture, we 
get 

Hi= £ ^^P(-iSt)[\+M+j\-\-i)(-i\ 

.7=1,2,3 

+ exp(i2fl j t)\+ j ){- j \ - cxp(-i2fl j t)\- :j )(+ : j\} + h.c, 

(6) 

where |±j) are the eigenstates of operator a Xj with eigenval- 
ues ±1. In the strong driving regime: fl S> 6, g, we can elimi- 
nate the fast-oscillating terms in Eq. (6) and then have l30ll32Tl 

Hi = y~] —<J Xj W ex P (-iSt) + a exp (iSt)] . (7) 

J=l,2,3 

Note that the operator set {a x . a x < , <fto Xi , aa Xj , 1} (j, j' = 
1,2,3, and j ^ j') form a closed Lie algebra, the time evolu- 
tion operator related to the above Hamiltonian can be formally 
written as 1I3TI1 

Ui(t) = exp[-iC{t)]Y[eaqp[-i(B s {t)aa Xj + B*{t)o) a Xj )] 
3 

x JJ exp[-iA jf (t)a x .a x .,}, (8) 
3^3' 



with the parameters determined by 

A jf (t) = ^f[^(exp(-iSt)-l)+t}, 
B j (t) = -^[e X p(iSt)-l], 

C(t) = J2%[^p(-iSt)-l) + t], (9) 

3 

and Ajj, (0) = Bj(0) = C(0) = 0. 

Suppose that all the qubit-cavity couplings are homoge- 
neous, i.e., gj = g (for j = 1,2,3) and set St = 2nir for 
integer n, we have B(t) = B*(t) = 0. Then, the time evolu- 
tion operator reduces to a simple form 

U I (t)=exp(-i^tS%), (10) 
with S x = Sj=i Vxj/2. Return to the Schroinger picture, 

U s (t) = U Q (t)Ui(t) 

= cxp (— iuia' at) J^J cxp (—iQ,j(T Xj t)Ui(t) 
3 

g 2 

= cxp (— iuxrat) cxp (— i2QS x t — i-^-tS x ).(ll) 

Here Oj = for gj = g mentioned above. Note that the 
effective coupling S x can be utilized to directly realize the 
multi-qubit GHZ state, once the relevant parameters are ap- 
propriately chosen ll33ll . Assume that the three-qubit register 
is initially at the state 

|V(0)) = |000), (12) 

where |1) (|0)) denotes the eigenstate of a z , <j z \1) = 1, 
c z |0) = — 1. Using the spin representation of atomic 
states for the operator S z , (S z = J2j=i u z 3 -/2)> the three- 
qubit states |000) and |111) can be expressed as collective 
states |3/2, —3/2) and |3/2, 3/2), respectively. Here, |J = 
3/2, M) is the eigenstate of the operators S z with the eigen- 
value M, M = — J, J. In terms of the eigenstates of 
S x El, we have 

3/2 

|3/2, -3/2) = £ c M \3/2,M) x , (13) 

M=-3/2 

and 

3/2 

|3/2,3/2)= J2 c M (-l) 3/2 - M \S/2,M) x , (14) 

M=-3/2 

where M = M'+l/2 and M' is an integer. As a consequence, 
the evolution of the system can be conveniently expressed as 
(up to a global phase factor) 

m)) = u s (t)\m) 

3/2 , 

2J c M exp [-i2fltM-i^j tM 2 } \3/2,M) x 

M=-3/2 

= -4(|3/2, -3/2) + i|3/2, 3/2)), (15) 
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with the choice g 2 t/5 = (4fc + l)n/2 and Vtt = (2TO + 3/4)vr 
(k, m are integers). Obviously, at t = T n the desired GHZ 
state JHQ] 

|V(T n )) = [/ S (T„)|^(0)) = -^=(|000) + *|111» (16) 

is obtained. In the above, the relations of the integers k, m 
and n are given by 

S 2 S 2 2S 2 3S 2 

n=^k + —r,n= m+- — . (17) 

9 4 S ^9 4eg 



B. Confirming the existence of the GHZ entangiement 

The GHZ state prepared above by one-step operation can be 
robustly confirmed by using the standard quantum-state tomo- 
graphic technique, i.e., reconstructing its density matrix. Such 
an approach was usually utilized to confirm the quantum state 
engineering in trapped-ions JUt, linear optics lfl7Ul9ll and the 
solid-state qubits I3a - l38ll . etc. However, these confirmations 
require many kinds of single-basis projective measurements 
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assisted by a series of quantum operations, and thus 2 
kinds of projections are needed for reconstructing a TV x TV- 
matrix, in principle. 

Fortunately, a significantly sim ple approach, i.e., spectral 
joint-measurements of the qubits 11391 14011 . can be utilized to 
high-effectively implement the desired confirmation. By this 
approach the states of the qubits could be jointly detected by 
probing the steady-state transmission spectra of the driven 
cavity, which is commonly coupled to the qubits. For the 
present case, the qubit-cavity detuning A = ui — u> r is as- 
sumed to be much larger than the coupling g (i.e., the system 
works in the dispersive regime) and the qubit-cavity couplings 
take the form H c = a^a^^ =1 Tj<j Zj . This indicates that the 
qubits cause the state-dependent frequency shift of the cavity. 
For example, if the qubits is prepared at the joint eigenstate 
1 000) (or | 111)) of the three qubits, then the frequency of the 

cavity is shifted as — J2j=i Tj (or Ylj=i ^j). Due to such a 
pull, the frequency of the transmitted photons through the cav- 
ity is shifted, which is dependent on the joint eigenstate of the 
qubits. Thus, the steady-state transmission spectra through 
the driven cavity can mark all the possible joint eigenstates 
of the qubits. Generally, unknown qubits should be denoted 
as a superposition of all the possible joint eigenstates of the 
qubits. As a consequence, the measured transmission spectra 
(aJa(u)d)}ss of the driven cavity may appear multiple peaks 
versus the driving frequency u>d (see the Appendix for the de- 
tailed derivation); each peak marks one of the possible joint 
eigenstates of the qubits, and its relative height corresponds 
to the probability of this state superposed in the three-qubit 
unknown state. 

Specifically, for the GHZ state prepared above the steady- 
state transmission spectra of the driven cavity should reveal a 
two-peak structure, see, e.g., Fig. 1(a) with the typical param- 
eters: (ri,r 2 ,r 3 ,K) = 2tt x (50,230,350, 1.69)MHz. To 
show clearly the simulated results for the testing, the parame- 
ters Tj could be adjusted by adiabatically tuning the qubit- 
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FIG. 1 : (Color online) Two spectral joint-measurements to confirm 
the GHZ entanglement: (a) directly for the GHZ state: (|000) + 
|lll>)/>/2, and (b) for the state (|000) + |011) + |101) + |110})/2 
generated after performing unitary operations on the GHZ state. 
Here, the parameters are selected as (Fi, F2, F3, k) = 2ir x 
(50,230,350, 1.69)MHz, and |0) = |000),|1) = |001),|2) = 
|010),|3) = |011),|4) = 1 100), 1 5) = |101),|6) = ] 110), and 
1 7) = | 111), respectively. 



transition frequencies. Desirably, the frequency-shift loca- 
tions: — Ti — T-2 — Ts, T\ + T2 + T3 and the relative heights of 
these two peaks: 0.5, 0.5, indicate that two joint eigenstates, 
|000) and 1 1 11), are superposed in the measured state with the 
same superposition probability 0.5. Of course, such a spectral 
joint-measurement result is just a necessary but not sufficient 
condition to assure the desired GHZ state, since a statistical 
mixture of these two joint eigenstates may also yield the same 
spectral distributions. To confirm the state |V>ghz) is indeed 
the coherent superposition of the states 1 000) and 1 1 1 1 ) , we 
need another spectral joint-measurement by using the quan- 
tum coherent effect. This can be achieved by first applying the 

unitary operation l\ 3 j=1 R V] (tt/4) = ]lj=i ex P ( i7rcr w/4) to 
each qubit, yielding the evolution 

I^GHZ) -> |V'GHz) = -R»i(^)-R» a (^)-R»s(^)I^GHz) 

= |(|000) + |011) + |101) + [HO)), (18) 

and then performing the spectral joint-measurement. It is ex- 
pected 09l l40ll that four peaks with the same relative height 
0.25 should be observed (see, e.g., Fig. 1(b)), if the prepared 
state is nothing but the desired tripartite GHZ state. However, 
if the prepared state is a mixture of the joint eignestates |000) 
and I 111), then eight peaks would be observed. 



III. TESTING TRIPARTITE MERMIN INEQUALITY BY 
SPECTRAL JOINT-MEASUREMENTS 

With the GHZ state, we now discuss how to test the tri- 
partite Mermin inequality (1) by jointly measuring the three 
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qubits, simultaneously coupled to a driven cavity. The test 
includes the following two steps. 

First, local parameters 6j(j = 1,2,3) are encoded into 
the generated GHZ state (16) by performing the single-qubit 
Hadamard-like operations 



|7>|6> |5>|4>|3>|2> |1>|0> 



|7>|6> |5>|4>|3)|2> |1>|0> 



Rj(0j) = R Z] {e ] /2)R X] ^/A)R Z] {-e j /2) 
J_ / 1 ie lB 
= x 2 \ " "' 1 



(19) 



Here, the typical single-qubit gates R Zj (9) = exp (ia Zj 9) and 
Rxj (9) = exp (ia Xj 9) can be relatively-easily implemented, 
see e.g. l29l l40ll . After these encoding operations, the gener- 
ated GHZ state I^ghz) is changed as 

= i[(l + e j(ei+e2+f,3) )|000) 
-<*s_ ie i(fli+9 a ))|ooi) 



le 

-i(9 2 +8 3 



ie 
— e 
— e 

i — ie 



i+«a))|010) 

m )\0ll) 
")\100) 



-i(«i+9 3 ) _ e i03 )|HO) 
l (e 1 +9 2 +e 3 )^| 111 ^_ 



)|101) 



(20) 



Second, we perform the joint projective-measurements to de- 
termine the required correlation functions E(6i, 62,63) for 
various combinations of these local variables. 

Experimentally, the above two steps can be repeated many 
times, and thus the correlation function can be determined by 



E(6 1: 6 2 ,6 3 ) = Pn+Poo + Poio + Pooi 
— -P011 — P101 — -Plio — Pooo- 



(21) 



Here, J2i, 3 ,k=o,i Rijk = 1 with Pyfc is the probability of 
the state IV'ghz) collapsing to the joint basis \ijk). With 
these projective measurements, various correlation functions 
required can be measured and then the tripartite Mermin in- 
equality (1) can be tested. Theoretically, the correlation func- 
tion can be easily calculated as 



P(#l, #2, 63) 



W ghz\Pt\iP" ghz) 
-cos(6»i +9 2 +63). 



(22) 



with the joint projective operator Pt = <r Zl <8 u Z2 ® 
a Z3 = I 111) <111| + |100)(100| + |010)(010| + |001)(001| - 
|011)(011| - 1 101) <101| - 1 110) <110| - |000)(000|. For the 
suitable choices of the local observables, e.g., {9i,9 2 , 63, 6[, 
6' 2 , 9' 3 } = {0, 7r/4, tt/2, 7r/4, 7r/4, 7r}, we have the ideal value 
of the Q -parameter in Eq. (1): 



Qi = V2 + 1 > 2. 



(23) 



This indicates that the inequality (1), namely 
Qi < 2, is violated. Furthermore, for the parameters 



I n it :A:A A 



(b), 



-800 -400 400 



*T It "rtiA" 11 



-400 400 800 



0.037 - 




fr"TJrii"lT 

-400 400 800 



(d) 



M HAM M. 



(» d - (» r (MHZ) 



-400 



400 800 

co (MHZ) 



FIG. 2: (Color online) Transmission spectra of the driven cavity 
versus the detuning for the evolved state |j/i"ghz) with the classi- 
cal variables {9 1 ,e 2 ,9s,0' 1 ,e' 2 ,9' 3 } = {0,7r/4,7r/2,7r/4,7r/4,7r}. 
Here, (a)-(d) correspond respectively to the parameters {6[, 62, #3}, 
{01,02,03}, {01,02,0s}, and {0i,0 2 ,0 3 }. With these spectral dis- 
tributions, the correlation functions required for testing the Mermin 
inequality (1) can be calculated. Other parameters of the system are 
the same as those used in Fig. 1. 



{01, 62,63, 6^62, 63} = {7r/4,0,0,37r/4, tt/2, tt/2}, the 
above tripartite Mermin inequality is maximally violated, i.e., 
Q l - 2y/2. 

Like in the usual tomographic reconstructions only one ba- 
sis, e.g., \ijk), is collapsed for one kind of projective measure- 
ment Pijk = \ijk) (ijk\. This implies that seven kinds of pro- 
jective measurements are required to complete the above joint 
projection Pt- However, by the spectral joint-measurements 
introduced in Refs. l39l 14011 . the probabilities Pijk(i,j, k = 
0, 1) can be determined simultaneously by just the spectral 
measurements of the transmission through the driven cavity; 
each peak of the transmission spectra marks one of the basis 
\ijk), and its relative height refers to the relevant probabil- 
ity Pijk- Specifically, for one set of classical variables {9\, 
62, 63, 6[, 6' 2 , 6' 3 } = {0, tt/4, tt/2, tt/4, tt/4, tt}, Figs. 2(a-d) 
show how the spectra of the driven cavity distribute (versus 
the qubit-driving detuning) for the state (20). For instance, 
four peaks, marking respectively the basis states 1 000) , |011), 
1 101), and 1 110) , are shown in Fig. 2(a). Their relative heights 
are equivalent: Pooo = Pni = Ploi = P110 = 0.25. Thus, 
the correlation function between three local variables can be 
easily calculated as 

£(•71-/4,0,0) = 1, 
£(tt/2,0,0) = 0.704, 
P(7r/4,7r/2,0) = 0.704, 
P(7t/2,7t/2,0) = 0. 

Consequently, the numerical experimental result of the Q- 
parameter is 



2.408 w V2 + 1 > 2, 



(24) 
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and thus the tripartite Mermin inequality is violated. Sim- 
ilarly, for another set of classical variables {0 l7 9 2 , 03,0 1; 
02>03> = W 4 ,0,0,37r/4,7r/2,7r/2}, Figs. 3(a-d) show 
all the probabilities of eight bases in the present three- 
qubit system. Again, the involved correlation functions 
are calculated as (E(6[, 9 2 , 63), E(6i,e' 2 ,e 3 ), E{6 1 ,9 2 ,9' 3 ), 
E(9[,9 2 ,9' 3 ))=(0.704, 0.704, 0.704, -0.704). As a conse- 
quence, 

Q e = 2.816 « 2V2 > 2, (25) 
and the Mermin inequality (1) is violated more strongly. 
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FIG. 3: (Color online) Transmission spectra of the driven cav- 
ity versus the qubit-drive detuning for the set of local variables 

{01,03,03,01,02,03} = {7r/4,0,0,37T/4,7r/2,7T/2}. Other pa- 
rameters are the same as those in Fig. 1. 



IV. DISCUSSION 

We have proposed a direct and experimentally-feasible 
scheme to test tripartite Mermin inequality with cavity- 
qubit system, wherein quantum state of three qubits with- 
out any direct interbit coupling is detected by measuring the 
dispersively-coupled cavity. We have numerically demon- 
strated that the local-variable-dependent probabilities of vari- 
ous bases superposed in the local-variable-encoded GHZ state 
can be directly read out by the cavity transmission. With these 
probabilities, various correlation functions on the local vari- 
ables of individual qubits are easily calculated, and conse- 
quently the violations of the three-particle Mermin inequal- 
ity are tested. Specifically, a few examples were utilized to 
numerically confirm the tests. Certainly, the present proposal 
could be generalized to test various Bell-type inequalities with 
more than three qubits in a straightforward way. 

Note that in our numerical-experiments little deviations 
exist between our estimated results and the ideal predic- 
tions. For example, if the local variables are set as {9i,9 2 , 
3 , 9[, 6' 2 ,9' 3 } = {0,tt/4, tt/2, tt/4, 71-/4, tt}, the values of Q- 
paramter given by our numerical experiments is Q e = 2.408, 



which deviates the ideal values Qi = %/2 + 1 with a quan- 
tity AQ = Qi - Q e = 0.006. Also, for local variables 
{9^62,83, e[,e' 2 ,9' 3 } = {tt/4, 0,0, 3tt/4, tt/2, tt/2} the in- 
equality (1) should be maximally violated with Qi = 2\/2, 
but our numerical experiment yields Q e = 2.816 = Qi — 
0.012. These deviations are due to the existence of the dissi- 
pation of the cavity, which yields various finite widths of the 
transmission spectra through the driven cavity. As a conse- 
quence, the relative heights of the measured peaks are lower 
than those of the ideal (5-type peaks. Therefore, the violations 
of the Mermin inequalities are less than those for the ideal 
cases. But, it is sufficient to show the violation of the Mermin 
inequality. 
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APPENDIX: TRANSMISSION OF A DRIVEN CAVITY 

In this appendix, the transmission spectrum of a three-qubit 
in a driven cavity is calculated in detail. The transition fre- 
quencies of the three qubits are denoted as u\, ui 2 and 013, 
respectively. We assume that the dispersive condition 



(1 . 9j_ 9j9j 



9j9r 



«1, j^j' = 1,2, 3, (Al) 



is satisfied for ensuring the effective dispersive coupling 
a Zj a) a between the jth qubit and the cavity. These conditions 
also ensure that the interbit interactions are negligible. Also, 
Aj = Lu r — ujj denotes the detuning between the jth qubit and 
the cavity, and Ajj> = luj — ujj the detuning between the jth 
and j'th qubits. 

In a framework rotating at u>d, the effective Hamiltonian of 
the qubit-cavity system is 



H 



(-6 + r 1( T 2l + T 2 a Z2 + r 3 o- Z3 )a t a 



UJ\ LU'2 UJ'£ ( * , 

+ y<r« + ycr 22 + y^3 + e ( a + a )> (A 2 ) 



L0 2 



0J3 



. 3 - gpAj, u,j = loj + Tj, (j = 1,2,3), and 5 



where Tj 

bJd — uj t is the detuning of the cavity from the driving. The 
master equation for the complete system reads 

q = -i[H,g} + K(aga^ -a)ag/2- go) a/2), (A3) 

where g is the density matrix of the qubit-cavity system. 

From the above master equation, the equations of motion 
for the mean values of various expectable operators are 



d(a t d) 
dt 



2elm(a) 



(A4a) 
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d(a) k . d(aa zl a Z2 cr Z3 ) K 

-jf = ( t6 -2^ a ' le J t = - -)(aa zl a Z2 a Z3 ) 

- iTt (a<7 Zl )- iT 2 (aa Z2 ) - iT 3 (aa Z3 ), (A4b) - ie(a zl a Z2 a Z3 ) - iT 1 (da Z2 a Z3 ) 



with 

d(aa Zl ) 



iT 2 (aa Zl a Z3 ) - iT 3 {aa Zl a Z2 ), 

(A4i) 



(iS - -)(aa Zl ) - ie(a Zl ) - iT 2 (dcr Zl cr Z2 ) 



dt y 2' 

- iTi(a) -iT 3 (aa Zl a Z3 ), (A4c) 

r* v •/ v -r/- \ d(^^d{a^ = d{a^ = 

j = {t8--){aa Z2 )-te{a Z2 )-tT x {aa Z2 a Zl ) dt dt dt J ' 

- iT 2 (a) -iT 3 (aa Z2 a Zs ), (A4d) 



d(c 



dt 



(i6--)(aa Z3 )-ie(a Z3 )-iT 1 (aa Z3 a Zl ) d{o zl o Z2 ) = d(a zl a Z3 ) = d{a Z2 a Z3 ) 

iT 3 (a) — iT 2 (aa Z3 a Z2 ), (A4e) dt dt dt 



0, (A4k) 



and 



d(aa Zl a Z2 ) k . 

- [id - -){aa zl a Z2 ) - ie(a Zl a Z2 ) 



dt v 2 

— iTi (acr 3 

- iT 3 (aa Zl a Z2 a Z3 ), (A4f) 



d(<Tzi l Z2g23) = 0, (A41) 
dt 



dlafj a \ K Th e steady-state distribution of the intracavity photon num- 

Z2 Zs = (iS — -z){da Z2 a Z3 ) — ie(a Z2 a Z3 ) ber can be obtained by solving the Eqs. (A4 a-i) under the 

steady-state condition, i.e., all the derivatives in the left sides 
* 2\a&z 3 ) i 3\ a c Z2 ; of above equations equate 0. Then, by numerical method, the 

— iri(aer Zl <72 2 <7 Z3 ), (A4g) steady-state average photons number inside the cavity can be 

obtained. Similar to the single-qubit and two-qubit cases in 
.„ . Ref. J39ll40ll , information of these eight basis states in arbi- 

\aa Zl a Z3 j _ ^ —){acr zl (Tza) ~ ie(cr Zl cr Z3 ) trary three-qubit state can be extracted from the spectra of the 

cavity transmission, since each peak marks one of the eight 
iTi{do Z3 ) — iT 3 (da Zl ) bases, and its relative height refers to its probability super- 

iT 2 (aa Zl a Z2 a Z3 ) , (A4h) posed in the measured three-qubit state. 



dt v 2 
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